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Abstract. Decaying three-dimensional (3D) turbulence is studied via direct 
numerical simulations (DNS) for an isotropic non-rotating flow and for rotating flows 
with and without hclicity. We analyze the cases of moderate Rossby number and large 
Reynolds number focusing on the behavior of the energy spectrum at large scales and 
studying its effect on the time evolution of the energy and integral scales for E(k) ~ fc 4 
initial conditions. In the non-rotating case we observe the classical energy decay rate 
£-io/7 an( j a g row t n f the integral length proportional to t 2 / 7 in agreement with the 
prediction obtained assuming conservation of the Loitsyanski integral. In the presence 
of rotation we observe a decoupling in the decay of the modes perpendicular to the 
rotation axis from the remaining 3D modes. These slow modes show a behavior similar 
to that found in two-dimensional (2D) turbulence whereas the 3D modes decay as in the 
isotropic case. We phenomenologically explain the decay considering integral conserved 
quantities that depend on the large scale anisotropic spectrum. The decoupling of 
modes is also observed for a flow with a net amount of helicity. In this case, the 3D 
modes decay as an isotropic fluid with a constant, constrained integral length, and the 
2D modes decay as a constrained rotating fluid with maximum helicity. 
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1. Introduction 

Controversy on the invariance of integral quantities in decaying turbulence has aroused 
during the last years pQ. Integral invariants are required to build phenomenological 
theories to explain the decay of kinetic energy and the evolution of integral length scales. 
For isotropic and homogeneous turbulence, the conservation of the Loitsyanski integral 
/ for an initial spectrum E(k — > 0) ~ Jfc 4 was called to derive the classical energy 
decay rate E ~ t~ 10 / 7 and a growth of the integral length L proportional to t 2 / 7 [2\. In 
a similar fashion, for an initial spectrum E(k — > 0) ~ Sk 2 , the assumed conservation 
of the integral quantity S associated to the conservation of linear momentum leads 
to E ~ t~ 6//5 [3j. In practice, theses quantities were shown to be only approximately 
conserved in closures [I] and in numerical simulations [5, 6j depending on the large scale 
spectrum of the initial conditions. 

Less is known about the decay of turbulent flows in rotating reference frames. 
These flows have been largely studied due to its relevance for scientific and engineering 
problems. Their importance has motivated numerous theoretical, experimental, and 
numerical works. Applications are broad, including areas as diverse as turbo machinery 
and rotor-craft, convective region of the sun and stars, large-scale flows in oceans, and 
convective scales in the atmosphere. 

It is well known that solid-body rotation inhibits the non-linear direct cascade of 
energy toward small scales reducing the dissipation rate of kinetic energy in comparison 
with non-rotating flows. This reduced dissipation has been observed in simulations 
El, experiments P HO], and studied theoretically 111]. An increase in the integral 
length parallel to the rotation axis with time has also been reported [T2l [T3] for these 
flows. 

Resonant wave theory has been used to take into account the effect of rapid rotation 
in turbulence [TU [151 EE Ej • According to this approach, the energy is transferred from 
small to large scales by resonant triadic interactions of inertial waves. The theory also 
argues that the resonant interactions are responsible for driving the flow to a quasi-two- 
dimensional state. In resemblance with the classic Taylor- Proudman theorem for steady 
flows [14] , this result is often called the "Dynamic Taylor-Proudman Theorem" (see eg. 
[18]). and leads to the decoupling of slow modes which behave as an autonomous system 
of two-dimensional (2D) modes for the horizontal velocity components (perpendicular 
to the rotation axis) for strong rotation [T6l [T9] . 

As a result of this reduced non-linear coupling and dissipation, for decaying flows 
in the laboratory different scaling laws were observed as the flow decays [HI [20], from 
classical non-rotating values at small times changing to different power laws after a time 
of the order of In [21J, for example, an initial isotropic decay is observed, followed 
by a cross-over for R Q « 0.25 after which the energy decays slower (E(t) ~ t -3 / 5 ). Such a 
decay law was proposed by [8] based on the assumption of energy tranfer being governed 
by the linear time Q" 1 . Strong correlation of the vertical flow leading to the growth 
and subsequent saturation of the integral length by vertical confinement was observed 
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in [TO! [T21 ET] • This saturation was observed at a time proportional to Q~ 5 ^ 7 in [10J. In 
[22] large-scale columnar-structure formation through linear inertial wave propagation 
was observed. Large scales form columnar eddies aligned with the rotation axis and a 
linear growth of the axial integral length takes place once the Rossby number passes 
below a certain threshold (R a ~ 0.4). With the increase of rotation, energy is retained 
by stable large scale structures and prevented form cascading to small scales. In some 
cases, energy was observed to decay faster for larger rotation frecuency. 

In simulations, [231 EI] reported depletion of the nonlinear energy cascade and 
growth of anisotropy. Also, an increase on the energy decay rate with rotation 
frequency was observed for the isotropic as well as the perpendicular modes. Two- 
dimensionalization was reported in several works (see e.g., [25]), together with the 
formation of columnar structures [26] as seen in experiments. In |27J, three distinct 
regimes were observed depending on the rotation frecuency. At low rotation rates the 
flow behaves as non-rotating. At intermediate rotation rates, a strong coupling between 
rotation and non-linear interactions dominates (with a slower decay of the energy), 
and at high rotation rates viscous effects are dominant, damping the nonlinear effects. 
Recently, the cases of helical and non-helical rotating decaying turbulence with the 
integral scale of the size of the box were numerically studied in [28] , where it was found 
that the presence of net helicity decreases even further the decay rate of energy (see 
also [29J). Rotating helicoidal flows have applications in atmospheric research, helical 
convective storms being an example [30] . 

To explain some of the experimental and numerical results, an extension to 
phenomenological predictions for rotating turbulence for E(k — * 0) ~ Sk 2 and 
E(k —>■ 0) ~ Ik 4 initial spectra (usually known as Saffman and Batchelor spectra 
respectively) has been proposed [8]. It includes a slow-down factor of the energy flux 
due to the presence of Rossby waves involving two different timescales: a long timescale 
representative of the turbulence evolution and a short one associated with the rotation 
frequency Q. Conservation of S and I is then called to derive the asymptotic decay of 
energy for both initial spectra resulting in E ~ t~ 3 / 5 and E ~ t~ 5 ^ 7 respectively. In 
the case of constrained turbulence, phenomenology leads to E ~ t~ l . However, these 
phenomenological arguments do not consider the effect of anisotropies in the integral 
conserved quantities or in the decay laws. 

In this work we numerically study the decay of rotating helical and non-helical 
turbulence with an emphasis on the anisotropies that arise when rotation is present, and 
on how integral quantities may be modified. The paper is divided in five sections. In 
section 2 we introduce the equations and describe how we solved them, with information 
regarding initial conditions. In section 3 we consider as an example a non-rotating flow, 
which behaves in agreement with previous results. In sections 4 and 5 we show and 
analyze results for the rotating non-helical and helical cases respectively. In the presence 
of rotation, we observe a decoupling of the energy decay rates for the 2D and three- 
dimensional (3D) modes. Studying the low wave number behavior we propose that the 
conservation of two-dimensional integral moments may explain these decays. In section 
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6 we finally summarize the results. 
2. Numerical Simulations 

The Navier-Stokes equation for an incompressible fluid in a rotating frame is solved 
numerically. When rotation is present, the equation reads 

d t u + wxu + 2Cxu = - VP + z/V 2 u, (1) 

together with the incompressibility condition 

V • u = 0. (2) 

Here u is the velocity field, uj = V x u is the vorticity, V = p/p — (fl x r) 2 /2 + u 2 /2 
is the total modified pressure, p is the (unit) density, and v is the kinematic viscosity. 
We chose the rotation axis in the z direction so that ft = Qz, Q being the rotation 
frequency. Our integration domain is a cubic box of length 2tt with periodic boundary 
conditions and the equations were solved using a pseudo-spectral method with the 2/3- 
rule for de-aliasing. All runs were performed with a resolution of 512 3 grid points. The 
initial conditions were a superposition of Fourier modes with random phases and an 
energy spectrum E(k) ~ k A with modes randomly distributed in a spherical shell of 
wave numbers between 1 and 14. Details of the simulations are given in Table 1. 

Table 1. Parameters used in the simulations: t* refers to the time of maximum enstrophy; 
Re, R , R™, H and h are respectively the Reynolds, Rossby, micro-Rossby numbers, the total 
helicity and the relative helicity. All quantities are calculated at time t*. A resolution of 512 3 
grid points was used in all runs. 



Run 


V 


n 


Re 


R 


K 


H 


h 


t* 


1 


8.5 x 1(T 4 





420 






0.01 


1 x 10" 4 


0.6 


2 


8.5 x 1(T 4 


10 


450 


0.1 


0.95 


0.05 


4 x 10~ 3 


0.7 


3 


8.0 x 1(T 4 


10 


530 


0.07 


0.7 


6.5 


0.5 


1.5 



In table 1, Re = LU/u is the Reynolds number, = U/(2LQ) is the Rossby 
number based on the integral scale L, R% = w/(2Q) is the micro-Rossby number with 
to = (a; 2 ) 1 / 2 , H — (u • u>) is the total helicity, h = H/(\u\\u)\) is the relative helicity, 
and t* is the time when maximum enstrophy is attained and when theses quantities are 
measured. We use L defined as 

L = 2nE- 1 J E{k)k~ l dk. (3) 

Note that R™ and R% are one order of magnitude apart. This is required for rotation 
not to completely damp the non-linear term in the Navier-Stokes equation leading to a 
pure exponential decay (see [17] for more details). 
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Figure 1. Energy decay for the non- rotating case (run 1). After a transient, the 
self-similar decay agrees with classical Kolmogorov theory. The t~ 10 ^ 7 slope is shown 
as a reference. 



3. Non-rotating flow 

3.1. Phenomenological arguments 

The classical Kolmogorov phenomenology leads to the well known energy spectrum 

E(k) ~ e 2 / 3 AT 5 / 3 , (4) 

which for a decaying self-similar flow with E(t) ~ kE(k) and using the balance equation 
dE/dt ~ e, gives the result 

dE E 3/2 

i ~ — < 5 > 

L can depend on time and extra hypothesis are required to obtain the energy decay. 
If L ~ L (where L is the size of the simulation domain), then dE/dt ~ E 3 ^ 2 /L 
and it follows that 

E(t) ~ t~ 2 . (6) 

If L ^ Lq but the spectrum at large scales is ~ k A , conservation of an integral 
quantity can be assumed to derive the decay rate of this flow. Traditionally the 
conservation of an initial k A dependence in the low wave number spectrum E(k) has 
been related to the invariance of the Loitsyanski integral / which we define as: 

1= / r 4 (u-u>, (7) 
Jo 
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where (u • u') is the isotropic two-point longitudinal correlation function which depends 
solely on r. If conserved, from dimensional analysis it follows that / ~ L 5 U 2 , then 
dE/dt ~ E 17 / 10 // 1 / 5 and we finally get 

E{t) ~ r 10 / 7 (8) 

as obtained by Kolmogorov [2]. 

In practice, / evolves slowly in time and is only approximately conserved for an 
initial large scale spectrum ~ k 4 . If the large scale spectrum is ~ k 2 , then another 
integral quantity is approximately conserved [3], which leads to a decay E[t) ~ t _6//5 . 
In the next section we present a simulation (run 1 of table 1) which approximately 
follows the Kolmogorov decay (see also [5]). The rotating cases in section IV and V have 
the same large scale energy spectrum. Conditions where the integral length saturates 
(reaching the box size in numerical simulations) have been reported in [28J. 



3.2. Numerical results 

In run 1 the energy spectrum (not shown) peaks initially at k — 14 and maintains an 
approximately k 4 scaling for low wave-numbers. The time history E(t) for run 1 is 
shown in figure [TJ After an initial transient of about six turn-over times, it shows a 
self-similar decay that is consistent with the t~ 10 ^ 7 law. 

In order to test further the behavior at large scales, we calculate the Loitsyansky 
integral I for this isotropic flow. In a recent work Ishida et al. estimated / by fitting 
E = Ik 4 /24:7c 2 to the energy spectrum at large scales [5]. In their simulations (with 
spatial resolution of 1024 grid points and an initial peak of the spectrum near k = 40 
or k = 80) the interval where E ~ k 4 holds is large enough for them to do the fitting. 
In our case (512 grid points) this interval is shorter and the fitting is not possible. 
Consequently, we checked spectral isotropy and then estimated I using equation (j7]). 
The two-point longitudinal correlation function can be estimated in the isotropic case 
using p]: 

roc 

(u.u')(r) = 2 / E(k)(sinkr — krcoskr) / (kr) 3 dk, (9) 
Jo 

where E(k) is the isotropic energy spectrum. 

The evolution of I(t) / 1(0) is shown in figureEl Although I(t) decays monotonically, 
after a transient its evolution is slow and it decreases only to approximately half its 
maximum value after 60 turn-over times. In [5] it was shown that its conservation 
improves as the extent of the large scale ~ k 4 spectrum is increased. 

The approximate invariance of / can further be used to estimate the growth of L. 
Writing L ~ J 1 / 5 ^ 2 / 5 an replacing in equation (jBJ) we get 

L ~ t 2 l\ (10) 

Figure [2] shows the evolution of the integral scale L z based on the one-dimensional 
spectrum for run 1 calculated as: 

L z = 2-kE- 1 [ E(k z )k~ x dk. (11) 
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Figure 2. One-dimensional integral length scale L z for run 1. After a transient period 
of time L z grows approximately as t 2 ^ 7 as derived phenomenologically. 




After an initial transient L z asymptotically settles down in the simulation to a growth 
rate close to t 2 ^ 7 . 

The results presented so far are consistent with Kolmogorov theory for decaying 
homogeneous and isotropic turbulence where initial conditions allow for the integral 
length to grow. In the next section we consider the analogy for the rotating case. 
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Figure 4. Energy as a function of time for run 2. After a transient, a decay close to 
although shallower than t~ 5 ^ 7 is attained. A t~ 5 / 7 slope is plotted for reference. 



4. Rotating flow 

4-1. Phenomenological arguments 

In this section we analyze a flow subjected to solid-body rotation in the z axis with 
rotation frequency Q (run 2). In non- helical rotating turbulence a 



E(k) ~ e 1 ^ 1 ^-* 



(12) 



spectrum is typically assumed [Ml ED [32J, [331 El]- From the balance equation this 
spectrum results in 

dE/dt ~ {E/Lfl/n. (13) 

Again, there are at least three possible scenarios. For L ~ Lq, dE/dt ~ E 2 [28] and 

E(t) ~ r 1 . (14) 

When L ^ Lq and ^(fc) ~ A; 4 at large scales as in our runs, constancy of / ~ U 2 L 5 
can be used again so that dE/dt ~ E 12 ^ 5 / (I 2 ^ 5 Q) and [8] 

£(t) ~ r 5/7 . (15) 

Invariance of J also leads to 

L~t 1/7 . (16) 
Finally, details of the rotating case with E{k) ~ /c 2 can be found in [S]. 
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Figure 5. Energy spectrum E(k,t) for run 2 from t = 1 to t = 45 in steps of At = 1. 
Energy piles up at large scales shallowing the spectra. A k 4 slope is shown for reference 
only. 

4-2. Numerical results 

Figure H] shows the evolution of E(t). After an initial nearly inviscid period, a transient 
period leads to a decay rate slightly steeper than E ~ t~ 5 / 7 . 

The conservation of / assumed to derive E{t) ~ t~ 5 / 7 is associated to a preserved 
A; 4 spectrum at large scales. However, in the presence of rotation an inverse cascade of 
energy develops shallowing the spectrum at large scales as the energy piles up at low 
wave-numbers. This effect is visible in figure [5] where the isotropic energy spectra is 
plotted for different times. As a consequence / is no longer approximately conserved 
and varies fast. Figure [6] shows the evolution of I(t)/I(0) departing from a constant 
value and increasing by an order of magnitude in clear contradiction with the hypothesis 
of constant /. 

When rotation is present, isotropy breaks down, the flow becomes axisymmetric 
and a privileged direction (the z axis) exists. In figure [7] we show the evolution of the 
energy in the slow 2D modes of the velocity E(k\\ = 0,t) (where kn = k z are the wave 
numbers in the direction parallel to the rotation axis), together with the energy of the 
remaining 3D modes with E(k\\ ^ 0,t). The k\\ ^ modes dominate initially having 
an amplitude one order of magnitude larger than the modes with kn = 0. This initial 
dominance is a result of the choice of initial conditions where Fourier modes are excited 
randomly all over shells in Fourier space, which results in most of the energy in modes 
with fcii 7^ 0. 

As time evolves, energy from 3D modes is transferred to the kn = plane and, after 
approximately six turn-over times, there is a cross-over after which energy in the 2D 
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Figure 6. Isotropic Loitsyansky integral / for run 2 
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Figure 7. E(t) for 3D modes with fc|| ^ (solid) and for 2D modes in the fc|| = 
plane (dashed) for run 2. Note an apparent decoupling in the decay rate between the 
energy in the 2D and the 3D modes. 
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modes prevails. Thereafter, 3D modes decay as in an isotropic flow without rotation 
following E(kn 7^ 0) ~ t~ 10 / 7 while perpendicular modes follow a shallower decay close 
to but now shallower than E±_ ~ t -5 / 7 . This suggests that after t ~ 6 the evolution of 
the 2D and 3D modes decouple and they decay separately with their own decay rates. 
The Rossby Number at that time is R D ~ 0.015, and the behavior is consistent with 
predictions of wave turbulence theory that obtains a decoupling for very small Rossby 
number [191 EH] with the 2D evolution of the modes described by the 2D Navier-Stokes 
equation. 



Jf.,3. Phenomenology revisited 

This behavior naturally leads us to review two-dimensional integral moments in 
addition to the isotropic Loitsyansky integral already introduced. For two-dimensional 
turbulence, [6] and [35] suggest that three canonical cases exist: E(k — > 0) ~ Jk^ 1 , 
E(k — > 0) ~ Kk and Eik — > 0) ~ hok 3 , where J, K and I^d will be respectively 
defined here as 

J = J (w.w')rdr, (17) 
K = J (u.u')rdr, (18) 

and 



'2D 



J r 3 (u.u')dr. (19) 

Moreover, J and K are integral invariants of motion with conservation of K being 
associated with conservation of linear momentum and invariance of J a consequence 
of vorticity conservation. In the 3D rotating case, since the k\\ = modes are the 
ones approximately decoupled for R Q 1, and the equations for these modes are 
equivalent to the 2D Navier-Stokes equations [191 EH]) we ma y wonder whether these 
integral quantities are conserved (or at least evolve slowly with time) in that manifold. 
In the rotating case, the relevant increments are then r = r±_ with perpendicular to 
f2, and the associated wave vectors are k^. 

We start by showing in figure [H] the energy spectrum for the vertically-averaged 
velocity field, that is to say, for wave numbers k\ = k^ + ky (hereafter, the perpendicular 
spectrum). This spectrum maintains a form proportional to k\ (albeit slightly 
shallower) . 

In order to find slowly varying 2D-like integral quantities in the simulation, we 
calculate the time evolution of K and I^d inverting the following equation which follows 
from assuming axisymmetry PQ: 

E(k ± ) = J ^(u.u')k ± r ± J Q (k ± r ± )dr, (20) 

so that 

(u.u')(r ± ) = J 2E{k ± )J {k 1 _r ± )dk ± (21) 
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Figure 8. Evolution of the two-dimensional perpendicular spectrum E(k±,t) for run 
2 from t = 0.5 to t — 45 in steps of At = 0.5. Note a slightly shallower than fc 3 scaling 
for low wave numbers. 
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Figure 9. Normalized K (solid) and I^d (dashed) as a function of time for run 2. Both 
parameters have a slow, almost constant behavior in time, showing little variation. 



can be estimated from the perpendicular spectrum. The results for K(t)/K(0) and 
^2d(^)/-^2d(0) are plotted in figure Both magnitudes behave in a similar fashion, 
showing slow variations with a relative constant value over the simulated time. 

Invariance of K or I 2 d leads to different energy decay rates. For constant K we 
can write K ~ L\U]_Lo\\ (where Lqii is the size of the box in the direction parallel to 
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Figure 10. Energy decay for perpendicular (fen = 0) modes compensated by t~ a for 
a = 5/7 (solid), a = 2/3 (dot-dashed) and a = 1/2 (dashed); a = 2/3 adjusts our 
data better. 



O) and assuming the slow-down factor in the energy dissipation rate by waves, as done 
in equation (iTBl . dE^/dt ~ (E±/ L^) 2 /O. Replacing L±_ in the last equation we get 
dE ± /dt ~ E\L Q \\l{Ktt) leading to 

E ± ~ t- 1 / 2 . (22) 

For constancy of I2D we have I2D ~ £jL^i.A)|]> an d using the same arguments, from 
dEjdt ~ (E ± /L ± ) 2 /n we get dE ± /dt ~ E 5 / 2 L 1 ^ / (1^ fi), which finally leads to 

£jl ~ t~ 2 / 3 . (23) 

In order to see whether any of these decay laws adjust our data better than the 
isotropic ~ t~ 5 ^ 7 law, we plot the 2D energy evolution compensated by i _5//7 , t _1,/2 , and 
r 2 / 3 (figure CUD. The decay law E ± ~ t" 2/3 seems to adjust better our simulation, which 
is consistent with the slow variation of I2D and with the initial perpendicular spectrum 
close to k\. 

The integral length parallel to the rotation axis L\\ = L z also behaves as in the 
isotropic case. In figure [EH we plot the evolution of L\\ compensated by the laws for 
the isotropic and axisymmetric cases already deduced. Clearly, the isotropic t 2 / 7 scaling 
adjusts better the results, in agreement with the isotropic-like decay of the energy in 
the 3D fast modes. 

5. Helical Rotating Flow 

We finally discuss briefly the effect of helicity upon the decay rate of energy in a rotating 
fluid. In order to incorporate net helicity into the flow we use a superposition of Arnold- 
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Figure 11. Integral length parallel to the rotation axis for run 2 compensated by t~ a 
for a = 2/7 (solid) and a = 1/7 (dashed); t 2 ^ 7 gives a better agreement to our data. 

Beltrami-Childress (ABC) [36] initial conditions achieving an initial relative helicity 
h ~ 0.99. As in run 2, the ABC flows were added in all shells in Fourier space between 
wave numbers k = 1 to 14 with an isotropic spectrum ~ k A . 

Figure [12] shows a comparison between the energy as a function of time for the 
helical and non-helical rotating flows. The helical energy decay, shown in dashed line, is 
slower than the non-helical case. This retard has been associated with an inhibition of 
the non-linear transfer of the energy toward smaller scales by a direct cascade of helicity 
|40j . The behavior has also been observed in [28], where rotating helical and non-helical 
flows with constant integral length were studied. 

The distinct evolution in the free decay of the helical flow can also be understood 
in terms of a phenomenological theory similar to those already presented. In this case, 
the direct transfer is dominated by the helicity cascade. Assuming maximal helicity we 
have [1Q] 

E{k±) ~ e 1/4 fi 5/4 ^ 5/2 . (24) 

Then it follows that 

dE±/dt ~ E^_/(L 6 ± Q 5 ). (25) 

It is unclear at this point whether in the phenomenological analysis we should 
separate the decay of the 2D modes from the 3D modes as in run 2, and whether the 
integral scale changes in time or not. Therefore, in figure [13] we plot the 2D and 3D 
energy decays for run 3 (same as figure but for the helical case). Note the anisotropic 
initial state with a relative excess of energy in the modes with kn = due to the ABC 
initial conditions used. After the first nearly inviscid transient, both sets of modes seem 
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Figure 12. Total isotropic energy decay for helical (dashed) and non- helical (solid) 
rotating turbulence. The results show a decreased decay rate in the presence of helicity. 
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Figure 13. Decay of the energy for 3D modes with k\\ ^ (solid) and for 2D modes 
in the k\\ = plane (dashed) for the helical rotating flow (run 3) showing different 
scaling laws. 
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Figure 14. Integral length scales parallel (solid) and perpendicular (dashed) to the 
rotation axis as a function of time. Both lengths saturate fast to an almost constant 
value near the simulation domain length (2tt). 

again to decouple and decay with different laws: ~ t~ 2 for the 3D modes and ~ t^ 1 ^ 3 
for the 2D modes. As in the non-helical case, the 3D modes decay faster, following the 
same law derived for an isotropic non-rotating flow where L ~ L . Indeed, in this run 
the integral scales grow fast during the transient and reach lengths close to the size of 
the box L (= 27r) before the self-similar decay starts. This is illustrated in figure [HJ 
that shows the evolution of the parallel and perpendicular integral scales. After t ~ 6 
L\\ is almost saturated, and L±_ keeps growing slowly but close to its maximum value. 
This results from a fast inverse transfer of energy (see the evolution of the isotropic 
energy spectrum in figure [15]) that may be associated to the large amount of energy in 
the ku = modes in the initial conditions. 

The fast increase and saturation of L\\ and L± give as a result the decay of the 2D 
and 3D modes as in constrained turbulence. For the 3D modes the ~ t~ 2 decay then 
follows. For the 2D modes, using equation (25) and the approximate constancy of the 
integral lengths, we get a decay 

E x ~ r 1/3 (26) 

in agreement with the simulation. The study of the cases where the integral scales are 
not constant are left for future work, and may require the identification of anisotropic 
integral conserved quantities as in the previous section. 
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Figure 15. Isotropic energy spectrum E(k,t) for run 3 from t = 1 to t = 45 in steps 
of At = 1. 




Figure 16. Perpendicular energy spectrum E(k±,t) for run 3 from t = 1 to t = 45 in 
steps of At = 1. 
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6. Conclusion 

In this work we presented studies of decaying turbulence in the presence of rotation and 
helicity. A simulation (run 1) without rotation and helicity was used to introduce some of 
the phenomenological arguments in unconstrained decaying turbulence. The simulation, 
with a large scale spectrum close to ~ k 4 , shows a slowly varying Loitsyansky integral 
and a decay law consistent with Kolmogorov 's E(t) ~ t~ 10/7 law. Detailed studies of 
such a decay can be found in [5]. 

When extending these arguments to rotating turbulence, approximate conservation 
of isotropic integral moments (as e.g., the Loitsyansky integral) is often assumed 
(see e.g., [8]). A simulation of non-helical rotating turbulence (run 2) was shown to 
decay slightly faster than what is predicted by these arguments. We argued that the 
approximate decoupling of slow and fast modes predicted in wave turbulence theory for 
rotating flows at very small Rossby numbers leads to different decay laws for the energy 
in the 2D and 3D modes. The 3D modes decay in agreement with phenomenological 
predictions for isotropic and homogeneous turbulence, while the decay of the 2D modes 
is consistent with phenomenological results obtained assuming integral moments of the 
two-dimensional Navier-Stokes equation are approximately conserved. 

Finally, the effect of helicity in rotating turbulence was considered in run 3. Helicity 
decreases the decay rate of turbulence even further as the direct transfer is dominated 
by the direct helicity flux (see e.g., [lQl[28]), and helicity tends to decrease the amplitude 
of the non-linear term in the Navier-Stokes equation. The initial conditions considered 
led to the fast saturation of the integral scales, and as a result the 2D and 3D modes 
in this run decayed as constrained turbulence. The 3D modes decayed as in the non- 
rotating (constrained) case, while the 2D modes were observed to decay slower than 
what is predicted for constrained non-helical rotating turbulence and in agreement with 
predictions that consider the effect of helicity. 

The three simulations presented here are far from an exhaustive exploration of the 
possible decay laws that may develop in rotating turbulence, and a detailed study of the 
effect of changing the initial large-scale energy spectrum dependence is left for future 
work, as well as studies of the effect of initial anisotropies in the decay, and the effect 
of scale separation between the initial integral scale and the box size (see e.g., [5j for 
a study for isotropic and homogeneous turbulence), and parametric studies varying the 
Reynolds and the Rossby numbers. 
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